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In this talk I report on recent progress in a few areas closely related to the virtual 
Compton scattering studies. In particular, I discuss the quark-hadron duality 
estimate of the -y*p — » A+ transition, QCD sum rule calculation of the 77* — > 
7r° form factor, and application of perturbative QCD to deeply virtual Compton 
scattering at small t. 

1. Soft vs. hard dynamics in QCD 

QCD and virtual Compton scattering. The kinematics of the amplitude of 
the process j*p — > jp' can be specified by the initial nucleon momentum p, 
the momentum transfer r = p — p' and the momentum q of the initial virtual 
photon, q 2 = —Q 2 . The final photon momentum q' is then given by q' = q + r 
with q' 2 = 0. Other important momentum invariants are t = (p' — p) 2 = 
(q — q') 2 and s = (q + p) 2 . 

Taking Q 2 large, i.e., at least above lGeV 2 , one can hope to enter the 
region where the amplitude is dominated by short distances between the two 
photon vertices and pQCD may be applicable. In this situation, it is tempting 
to speak about the "virtual Compton scattering on a single quark" implying 
that the large-Q 2 behaviour is given just by the quark propagator (see Fig. la), 
while the long-distance information is accumulated in a distribution function 
F{X,t) described by the matrix element of (p'\q...q\p) type. However, the 
factorization formula 

M(Q 2 ,s,t) -> { m{s/Q 2 ,X)F(X,t)dX (1) 
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only makes sense if \t\ <C Q 2 - Otherwise, if \t\ ~ Q 2 , large momentum enters 
into the hadron wave function and one deals with the scattering process on the 
hadron as a whole. 




a) b) c) 



Figure 1: a) Handbag diagram contributing to the DVCS amplitude at t = 0. The lower 
blob corresponds to the double quark distribution F(x, y). b) One of the two-gluon-exchange 
diagrams dominating the VCS amplitude for asymptotically large Q 2 and t. c) Lowest- 
order diagram for the 7*7 — + tt° transition form factor. The blob corresponds to the pion 

distribution amplitude tp^ix). 



For \t\ < Q 2 , the function F(X,t) in eq.(pj) looks like a parton distribu- 
tion function f(x) with an additional form-factor-type dependence on t. To 
make analogy with deep inelastic scattering, it is instructive to recall that the 
imaginary part of the virtual forward Compton amplitude (for which q' = q 
and p' — p) in the limit of large Q 2 and fixed Bjorken variable ( = Q 2 /2(pq) 
can be written as 

f 1 f(x)6((q + xp) 2 ) ■ 2(qp)dx = f 1 f(x)5(x - Q 2 /2(qp))dx = /(C) , (2) 
Jo Jo 

where xp is the fraction of the initial hadron momentum carried by the in- 
teracting quark. The usual parton distribution functions /(C) correspond to 
exactly forward matrix elements, with r = p' — p = 0, while the kinematics 
of VCS requires that r and t = r 2 ^ 0. Hence, we need a new type of 
parton functions F(X, t)u. Inlhe limit t — > 0, they reduce to the "asymmetric 
distribution functions" F(X) on. Hence, the studies of deeply virtual Comp- 
ton scatterinelDVCS) are related to a new field of pQCD applications. As 
shown in refsErH, the asymmetric distributions F(X) have features of both the 
distribution functions and distribution amplitudes (wave functions). A more 
detailed discussion of DVCS at small t will be given in Section 3. 

Another situation in which pQCD is applicable is when both Q 2 and \t\ are 
asymptotically large. Then the virtual Compton scattering amplitude factor- 
izes into a convolution of the short-distance amplitude m({xi}, {j/j}, Q 2 , s, t) 
and two distribution amplitudes f {x\, X2, £3), <p{yi, 1/2,1/3) describing the pro- 
ton in the initial and final state, respectively (see Fig. 16). They are related to 



2 



matrix elements of (0\q . . . q . . . q\p) type and give the probability amplitude, 
that e.g., the initial proton can be treated as three collinear quarks with the 
momentum p divided into fractions xip, X2P, x^p with xx + X2 + £3 = 1. The 
short-distance amplitude m({xi}, {yj}, Q 2 , s, t) is given by Feynman diagrams 
involving two hard gluon exchanges, suppressed by a factor (a s /7r) 2 ~ 1/100 
compared to the "soft contribution" produced by a simple overlap of soft wave 
functions, without any gluon exchanges. The soft term, however, has an ex- 
tra power of 1/Q 2 for large Q 2 . As a result, the hard term asymptotically 
dominates, though the soft term may be much larger than the hard one for 
accessible Q 2 . 

Quark-hadron duality and j*p — > A + transition. In particular, a purely 
soft contribution to G P M (Q 2 ) calculated within the local quark-hadron duality 
approach □ is in good agreement with experimental data up to Q 2 ~ 20 GeV 2 . 
The same approach was used recently U to get the estimates of the soft term 
for the 7*p — > A + transition. For the magnetic form factor G* M (Q 2 ) these 
estimates are rather close to the results of the analysis of inclusive SLAC data 
□13. A small value for the ratio G* E (Q 2 ) / G* M (Q 2 ) obtained in rem is-also in 
agreement with available dataB, in contrast to the pQCD prediction □ which 
gives G* E (Q 2 ) / G* M {Q 2 ) — > —1 for the ratio of hard contributions. 
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Figure 2: Quark-hadron duality estimates for the 7*p - 



A+ transition. 



Left: form factor 



G* M (Q ). Right: ratio of form factors G*„(Q ) and G* M {Q ). 



Within different nonperturbative approachesEHlcfEilO, it was observed that 
soft terms are sufficiently large to describe the data or that the hard terms are 
too small compared to the data. Hence, there is a growing evidence that soft 
terms dominate the exclusive amplitudes at accessible energies. Of course, the 
magnitude of the hard contribution depends on the shape of distribution am- 
plitudes (DA's). The latter are usually integrated with the_syeights like l/xxX2, 
so the humpy DA's of Chernyak-Zhitnitsky (CZ) typeEjO produce contribu- 
tions which are much larger than those obtained with smooth DA's close to the 
"asymptotic" forms. The CZ wave functions were originally motivated by QCD 
sum rule analysisll3. However, the results of the QCD sum rule calculations of 
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the DA's are extremely mo del- dependent and unreliable. Furthermore, for the 
theoretically most clean case of the 7*7 — > ir° form factor, both a direct QCD 
sum rule calculation of this form factor and available experimental data show 
no enhancement compared to the pQCD result obtained with the asymptotic 
DA for the pion. 

2. 7*7 — > 7T° form factor 

pQCD analysis. The transition 7*7* — + tt° of two virtual photons 7* into a 
neutral pion provides an exceptional opportunity to test QCD predictions for 
exclusive processes. In the lowest order of perturbative QCD, its asymptotic 
behaviour is due to the subprocess 7*(g) + 7*(<?') — ► q(xp) + q(xp) with x (x) 
being the fraction of the pion momentum p carried by the quark produced 
at the q (q 1 ) photon vertex (see Fig.lc). The relevant diagram is similar to 
the handbag diagram for deep inelastic scattering, with the main difference 
that one should use the pion distribution amplitude ^(x) instead oflparton 
densities. For large Q 2 , the perturbative QCD prediction is given by til 

(3) 

Experimentally, the most important situation is when the lower virtuality pho- 
ton is (almost) real q' 2 f=a 0. In this case, necessary nonperturbative information 
is accumulated in the same integral I (see eq.([|)) that appeArsjn_the one-gluon- 
exchange diagram for the pion electromagnetic form factorEJ'EJO. The value of 
/ depends on the shape of the pion distribu tkua amplitude tp w (x). In particular, 
using the asymptotic form ^p^ix) = 6/ T xxt§liZl gives F™ t ^ (Q 2 ) = iirf^/Q 2 
for the asymptotic behaviourlij. If one takes the Chernyak-Zhitnitsky formlij 
ip^ z (x) — 30/ 7r x(l — x)(l — 2x) 2 7 the integral I increases by a sizable factor of 
5/3, and this difference can be used for experimental discrimination between 
the two forms. 

Note, that the pQCD hard scattering term for 77* — * ir° has the zeroth 
order in the QCD coupling constant a s , just like in deep inelastic scattering. 
Hence, there are good reasons to expect that pQCD for F 77 . 7r o(Q 2 ) may work 
at rather low Q 2 . The Q 2 = limit of F 11 ,-aiQ 2 ) is known from tt° — » 77 
decay rate. Using PCAC and ABJ anomaly 113, one can calculate Ry~*^o(0) 
theoretically: i 7 77 . 7r o(0) = l/nf n . It is natural to expect that a complete QCD 
result does not strongly deviate from a simple interpolation!!^ 7r/ 7r i 7 ' 77 , 7r o (Q 2 ) = 
1/(1 + Q 2 /4ir 2 f 2 ) between the Q 2 = value and the large-Q 2 asymptotics |J 

c In particular, such an interpolation agrees with the results of a constituent quark model 
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This interpolation implies the asymptotic form of the distribution papplitudc 
for the large^O 2 limit and agrees with CELLO experimental dataEJ. It was 
also claimedEj that the new CLEO data available up to 8 GeV 2 also agree with 
the interpolation formula. 

Comparing the data with theoretical predictions, one should take into 
account the one-loop pQCD radiative corrections to the hard scattering ampli- 
tude calculated in refE2l. Effectively, the correction decreases the leading-order 
result by about 20%, still leaving a sizable gap between the prediction based 
on the CZ amplitude-and the phenomenologically successful Brodsky-Lepage 
interpolation formulia Hence, the new preliminary datfa seem to indicate 
that the magnitude of I is close to that corresponding to the asymptotic form 
of the pion distribution amplitude. Because of the far-reaching consequences 
of this conclusion, it is desirable to have a direct calculation of the 77* — > ir° 
form factor in the intermediate region of moderately large momentum transfers 
Q 2 >1 GeV 2 . As we will see below, the QCD sum rules allow one to calculate 
F 77 . w o(Q 2 ) for large Q 2 without any assumptions about the shape of the pion 
distribution amplitude, and the result can be used to get information about 
ip„(x). 

QCD sum rules and pion distribution amplitude. The CZ sum rule written 
directly for the pion distribution amplitude v>tt(x) is 

1 8 7ro «^ . {ll[6{x) + 5(1 - x)\ + 2[5'(x) + 5'(1 - x)]}. (4) 
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Here, M is the auxiliary Borel parameter which must be taken in the region 
where the r.h.s. is least sensitive to its variations, and sq is the effective 
onset of the continuum_6tted to maximize the M 2 -stability region. From the 
QCD sum rule for /j^E-a sq w 0.7 GeV 2 . As emphasized in refc3, the lowest 
condensates (GG) and (qq) 2 taken into account in eq.(^) do not provide all the 
information necessary for a reliable determination of (p w (x). The humpy CZ 
shape is, in fact, a compromise between the S(x), 5(1 — x) condensate peaks 
and the smooth x(l — x) behaviour of the perturbative term. Adding higher 
condensates, e.g., (qD 2 q), one would get even higher derivatives of 5(x) and 
5(1 — x). The sum of such singular terms can be treated as an expansion of 
some finite- width function 5<p(x): 

8<p(x) = a 5(x) + a\8'{x) + a 2 5"(x) + . . . + {x 1 - x}. (5) 



calculation!; 
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Of course, the knowledge of eto alone is not sufficient for a reliable reconstruc- 
tion of 6tp(x). On the other hand, the higher coefficients a±, 02, etc. are given 
by a sum of several higher condensates whose magnitudes are completely un- 
known. Hence, no strict conclusions can be made. The CZ procedure is equiv- 
alent to assuming_that a\,a%, ... ~ 0, though other choices (e.g., nonlocal 



condensate modelEj) may look more realistic. 

QCD sum rule for doubly virtual form factor F~*y*„o(q , q' 2 ). Instead 
of following the steps dictated by the old logic: 1) pQCD factorization for 
i 7 7 , 77r o(Q 2 ); 2) QCD sum rules for the moments of f^(x) (which are unreli- 
able); 3) calculation of I = j} (p w (x)/xdx, we developed in ref. El the approach 
which 1) starts with the QCD sum rule for F 7 » 7 »^o (q 2 , q' 2 ) in the q' 2 — » limit; 
2) information about I is extracted from this sum rule and 3) then used to 
make conclusions about the shape of ip 7T (x). When both virtualities of the 
photons are large, we have the following QCD sum rule: 

nUF Tl ^(Q 2 ,q' 2 )^2 Tdse-^ 2 C , ^_ {xQ \ ' 7 ^ 



(6) 



Jo 
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In this situation, the pQCD approach is also expected to work. Indeed, ne- 
glecting the sxx-term. compared to xQ 2 — xq' 2 and keeping only the leading 
0(1/Q 2 ) and 0(1/ q' 2 ) terms in the condensates, we can write eq.(^) as 

, 9n 4tt f 1 dx f 3M 2 , , /M 2, 

f ^ q ) =w,L \^ (1 - e )xx 

+ ^(^GG)[S(x)+S(x)} 
,(qqf (ll[S(x) + S(x)} + 2[5'(x) + S'(x)]j } . (7) 



81M 4 

The expression in curly brackets coincides with the QCD sum rule (Q) for the 
pion distribution amplitude f„ip„(x). Hence, when both Q 2 and q' 2 are large, 
the QCD sum rule (||) exactly reproduces the pQCD result (ph. 

One may be tempted to get a QCD sum rule for the integral / by taking 
q' 2 = in eq.(|^). Such an attempt, however, fails immediately because of 
the power singularities 1/q' 2 , 1/q' 4 , etc. in the condensate terms. It is easy 
to see that these singularities are produced by the 8(x) and §'(x) terms in 
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eq.(Q). In fact, it is preciseb^these terms that generate the two- hump form for 
(Pir(x) in the CZ-approachtL3. The advantage of having a direct sum rule for 
F 77 « 7r o(Q 2 . q' 2 ) is that the small-o/ 2 behavior of i 7 77 » w o (Q 2 , q' 2 ) is determined 
by the position of the closest resonances in the q' channel, which is known. 
Eventually, l/q' 2 is substituted for small q' 2 by something like 1/m 2 and the 
QCD sum rule in the q' 2 limit is 



Q 2 - 2s 



2 Q 2 - 6s + 3s 2 /Q 2 ( s* p s p 
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In Fig.||, we present a curve for Q 2 K y7 * w o(Q 2 )/47r/ 7I . calculated from eq.(||) 
for standard values of the condensates, p- and 7r-meson duality intervals s p = 
1.5 GeV 2 <B, s = 0.7 GeV 2 and M 2 = 0.8 GeF 2 . It is rather close to the curve 
corresponding to the Brodsky-Lepage interpolation formula 7r/ w i 7 ' 77 « 7r o (Q 2 ) = 
l/(l + Q 2 /47r 2 / 2 ) and to that based on the p-pole approximation i:f v F{Q 2 ) = 
l/(l + Q 2 /m 2 ). Hence, our result favors a pion distribution amplitude which is 
close to the asymptotic form. It should be noted, that the p-pole behaviour in 
the Q 2 -channel has nothing to do with the explicit use of the p-contributions in 
our models for the correlators in the g' 2 -channel: the Q 2 -dependence of the p- 
pole type emerges due to the fact that the pion duality interval so ~ 0.7 GeV 2 
is numerically close to m 2 » 0.6 GeV 2 . Taking the lowest-order perturbative 
spectral density p quark (s, q a — 0, Q 2 ) = Q 2 /(s + Q 2 ) 2 and assuming the local 
quark-hadron duality, we obtain the result 



UF^oiQ 2 ) = - [ S ° p^ ark (s,0,Q 2 ) 
t Jo 



ds = 



tt(1 + Q 2 /s ) 



(9) 



coinciding, for so = 47r 2 / 2 w 0.67 GeV with the BL-interpolation formula. 
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Figure 3: Combination Q 2 i ? 77 » 7r o (Q 2 )/47r/ 7r as calculated from the QCD sum rule (solid 
line), p-pole model (short-dashed line) and Brodsky-Lepage interpolation (long-dashed line). 

Lessons. 1) CZ sum rule is an unreliable source of information about the 
pion DA; 2) Pion DA is narrow; 3) Since the diagrams for the nucleon DA's 
have the structure as in the pion case, we should expect that the nucleon 
DA's are also close to asymptotic and , hence, the two- gluon- exchange hard 
terms are very small for accessible Q 2 and t. Thus, it is very important to 
get the estimates for the soft contributions to the virtual Compton scattering 
amplitude (see, e.g., refs. E3, where the high-i real Compton scattering on the 
pion was considered). 

3. Small-i, large-Q 2 limit of VCS: a new pQCD area 

Recently, X. Ji § suggested to use the deeply virtual Compton scattering 
(DVCS) to get information about some parton distribution functions inacces- 
sible in standard inclusive measurements. He also emphasized that the DVCS 
amplitude has a scaling behavior in the region of small t and fixed XBj which 
makes it a very interesting object on its own ground. 

Double distributions. In the scaling limit, the square of the proton mass 
m 2 = p 2 can be neglected compared to the virtuality Q 2 = —q 2 of the initial 
photon and the energy invariant p-q = m v v. Thus, we set p 2 = and, for small 
t, we also have r 2 = 0. Then the requirement p' 2 = (p — r) 2 = p 2 reduces to the 
condition p ■ r — which can be satisfied only if r is proportional to p: r — (p, 
where C coincides with the Bjorken variable XBj = Q 2 /2{p • q)i < xbj < 1- 
Naturally, the light-like limit of 4-momenta p, r is more convenient to visualize 
in a frame where the initial proton is moving fast, rather than in its rest frame. 

Though the momenta p and r are proportional to each other, one should 
make a clear distinction between them since p and r specify the momentum 
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flow in two different channels. Since the initial quark momentum originates 
both from p and r, we write it as xp + yr. In more formal terms, the relevant 
light-cone matrix elements are parameterized as 

(p-r\ M0)zE(0, z;A)il> a (z)\p)\*=o = «(P ~ r)2u(p) (10) 
L So { e ~ iX{PZ) ~ iV{rZ)Fa ( x > y } ~ e ix(pz) ~ mrz) Fu(%>y)) 0{x + y< l)dydx, 

etc., where z = 7 M z^ and u(p — r),u(p) are the Dirac spinors for the nucleon. 

Taking the limit r = gives the matrix clement defining the parton dis- 
tribution functions f a (x), fa{%)- This leads to the reduction formula: 

l—x 

F a (x,y)dy = f a (x). (11) 

Asymmetric distribution functions. Since r = (p, the variable y appears 
in eq.([l0|) only in x + yC, = X and x — yC, = X — £ combinations, where X and 
(X — £) are the total fractions of the initial hadron momentum p carried by 
the quarks. Integrating the double distribution F(X — y) over y we get the 
asymmetric distribution function 

,min{X/C,X/C} 

J?W= / F a (X-y<:,y)dy, (12) 

Jo 

where £ = 1 — C- Since C < 1 and x + y < 1, the variable X satisfies a 
natural constraint < X < 1. In the region X > ( (Fig. 4a), the initial 
quark momentum Xp is larger than the momentum transfer r = Cp, an d we 
can treat J-^{X) as a generalization of the usual distribution function f a (x). 
In this case, the quark goes out of the hadron with a positive fraction Xp of 
the original hadron momentum and then comes back into the hadron with a 
changed (but still positive) fraction (X — ()p. The Bjorken ratio £ specifies 
the momentum asymmetry of the matrix element. Hence, one deals now with 
a family of asymmetric distribution functions J-^(X) whose shape changes 
when C is changed. The basic distinction between the double distributions 
F{x,y) and the asymmetric distribution functions T^(X) is that the former 
do not depend on the momentum asymmetry parameter £, while the latter are 
explicitly labelled by it. When £ — > 0, the limiting curve for Tq{X) reproduces 
the usual distribution function: 

J? =Q (X) = f a (X). (13) 

Another region is X < £ (Fig. 46), in which the "returning" quark has a negative 
fraction (X — £) of the light-cone momentum p. Hence, it is more appropriate 
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Figure 4: Momentum flow corresponding to the asymmetric distribution function in two 
regimes: a) X > ( and b) X = Y( < (. 



to treat it as an antiquark going out of the hadron and propagating together 
with the original quark. Writing Iasl = Y(, we see that the quarks carry 
now positive fractions Y^p = Yr and Yr = (1 — Y)r of the momentum transfer 
r, and the asymmetric distribution function in the region X = Y( < £ looks 
like a distribution amplitude ^((Y) for a qq state with the total momentum 
r = (p: 

*cQ0 = / F((Y-y)C,y)dy. (14) 
Jo 

Leading-order contribution. Using the parameterization for the matrix 
elements given above, we get a parton-type representation for the handbag 
contribution at t = 0: 

T s % v mm (p, q, r) = (<T - -£7ZW + P V 1")] E e " a/W (?v(0 + ^(C)), 

a 

where only the {/i <-► ^-symmetric part is shown explicitly and Ty(£) is the 
invariant amplitude depending on the scaling variable (: 

™=£( Y ^ T - + ±)f<(X)dX. (15) 

The term containing 1/(X — ( + ie) generates the imaginary part: 

--ImT£ (0=^(0- (16) 

7T 

Though F"_ (X) = f a (X), in the general case when ( ^ 0, these two 
functions differ. Furthermore, the imaginary part appears for X — £, i.e., in a 
highly asymmetric configuration in which the second quark carries a vanishing 
fraction of the original hadron momentum, in contrast to the usual distribution 
/o(C) which corresponds to a symmetric configuration with the final quark 
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having the momentum equal to that of the initial one. A characteristic feature 
of the asymmetric distribution functions T^(X) is that they rapidly vary in 
the region X<(~ and vanish for X = 0. However, the limiting curve J-^=o(X) 
does not necessarily vanish for X = 0, i.e., the limits £ — > and X — > do not 
commute. For this reason, if £ is small, the substitution of FHX) by f a (X) 
may be a good approximation for all X-values except for the region X< £, and 
it is not clear a priori how close are the functions J 7 ? (£) and f a (C)- 

Evolution of the double distributions. The purely scaling behavior of the 
DVCS amplitude is violated by the logarithmic <5 2 -dependence o{Fns(%> V', Q 2 ) 
governed by the evolution equation 



d f f 

Q^ F Ns{x,y,Q 2 ) = d£ R N s {x,y; $,,77; g)F NS {£,,r); Q : 
dQ Jo Jo 



)dv (17) 



(the flavor-nonsinglet (NS) component is taken for simplicity). Since integra- 
tion over y converts F]\fg(x,y) into the parton distribution function /jvs(x), 
whose evolution is governed by the GLAPD equation E3'E3Ej, our kernel has 
the property 

l—x y 

RNs(x,y;£,n;g)dy = ~P NS (x/£). (18) 

For a similar reason, integrating Rns{x, y; £, rj; g) over x one should get the 
evolution kernel V(y, ij; g)\!3'tB for the pion distribution amplitude 

Rns(x, y; f , jj; g)dx = V(y, rj; g). (19) 

In the formal Q 2 — > oo limit, F(x,y;Q 2 — > oo) ~ 6(x)yy, i.e., in each of its 
variables x, y, the double distribution tends to the characteristic asymptotic 
form: 8(x) is specific for the distribution functions, while the yy-form is .±he 
asymptotic shape for the lowest-twist two-body distribution amplitudes E30. 

Evolution of asymmetric distribution functions. As a result, the evolution 
of the asymmetric distribution functions T?(X) proceeds in the following way. 
Due to the GLAP-type evolution, the momenta of the partons decrease and 
distributions become peaked in the regions of smaller and smaller X. However, 
when the parton momentum degrades to values smaller than the momentum 
transfer r = C,p, the further evolution is like that for a distribution amplitude: 
it tends to make the distribution symmetric with respect to the central point 
X = C/2 of the (0, C) segment. 

Conclusions. DVCS opens a new class of scaling phenomena characterized 
by absolutely new nonperturbative functions describing the structure of the 
proton. The continuous electron beam accelerators like TJNAF and ELFE 
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may be an ideal place to study DVCS 0. The asymmetric distributions can 
also be studied in the processes of large-Q 2 meson electroproductionH, etc. 
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